clc
clear
mpc=loadcase(casomasgrandebucle3)
n = [300 200 300;0 0 0;200 210 220;20 10 20;0 0 2;30 30 30;20 10 10;31 0 31;200 210 190]
for i=1:3
mpc.bus(:,3)=n(:,i);
mpc=loadcase(casomasgrandebucle3)
rundcopf(casomasgrandebucle3) % resuelve el opf.
define_constants;
resultado=rundcopf('casomasgrandebucle3');
Mult_Lagrange(:,i)=resultado.bus(:,LAM_P);
end
display(Mult_Lagrange)


